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Abstract. In asymmetric heavy ion collisions like dA or pA, particle produc- 
tion yields are different in the forward [d- or p-side) and backward (A-side) 
rapidity directions. The rapidity distribution reflects the geometry and phase- 
space distribution of nuclear matter. These properties may depend on the 
time evolution of the collision. Due to the smallness of the backward-forward 
differences, the rapidity asymmetry factor can be useful to quantify nuclear 
modification effects, like e.g. shadowing and the EMC effect. Our work is 
a survey of the nuclear modification factor and the rapidity asymmetries at 
RHIC energies. We analyze the rapidity dependence and the strength of the 
nuclear effects. We focus on the high transverse momentum region, and make 
predictions for the role of nuclear modifications and rapidity asymmetries for 
future experimental measurements at increasing absolute values of rapidity. 
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1. Introduction 



Measuring nuclear modifications is not only interesting by itself, but is crucial to 
the understanding of the hard-probe signature of the quark-gluon plasma (QGP) 
formed in high energy heavy ion collisions. The modifications can be determined 
experimentally by different methods in a wide kinematical region: (i) the nuclear 
modification factor, R'^j^,{px) can magnify the deviation caused by collective nu- 
clear effects relative to 'simple' nucleon-nucleon collisions; (ii) the (pseudo)rapidity 
asymmetry, YAsvmiPT) 1 measures differences of the hadron spectra between back- 
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ward and forward directions. The deviation from unity of the above two quantities 
originates in the geometrical and nuclear properties of the colliding system. 

Measuring Rj^j^, {pt) and Y^^ {pt) was an important task for past and present 
experimental collaborations [[H [^j and is going to remain at the center of attention 
of future experiments. In this paper we would like to emphasize the connection be- 
tween the two quantities and present recent results [IS HI [5] for upgraded detectors 
at RHIC and future experiments at the LHC. 

2. Quantifying Modifications in the Nuclear Medium 

The nuclear modification measures the effect of the collective nuclear forces in a pA 
or AA' collision, compared to an average nucleon-nucleon collision. The deviation 
can be measured as a ratio of the given hadron {h) spectra per nucleon in AA' to 
the NN collisions. The nuclear modification factor can be defined for any given 
pseudorapidity: 



^^^^^' " {Nun) EuA^'at/d^PT 
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Nuclear effects can make R'^^{pt,ti) grater or smaller than I, representing an en- 
hancement or suppression, respectively, relative to the NN hadron spectra. 

In asymmetric collisions, hadron production at forward rapidities may be dif- 
ferent from what is obtained at backward rapidities. It is thus of interest to study 
ratios of particle yields between a given pseudorapidity value and its negative in 
these collisions. The pseudorapidity asymmetry YAsym{PT) is defined for a hadron 
species h as 
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Let us consider the ratio of the backward and forward nuclear modification 
factors in dAu collisions for species h: 
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(3) 



where (Nbin) is the average number of binary collisions in the various impact- 
parameter bins, and it is given by the thickness function of the Glauber model. 
Therefore, {N^^ ) = {N^^ ). Furthermore, the pp rapidity distribution is sym- 
metric around y = 0. Thus, if the same backward and forward (pseudo)rapidity 
ranges are taken in both directions (i.e. |?7m,m| < |''?| < |??maa;|), then the pp yields 
cancel in eq. ([3]) and one obtains that the ratio defined in eq. ([3]) is identical to the 
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pseudorapidity asymmetry eq. ^: 

YAsynAPT) - R'^ipr) 



R'^AuiPT,V<0) 
RdAu(PT,V>0) 



(4) 



3. Nuclear Modification at Forward and Backward Rapidities 

We analyzed several types of shadowing parameterizations such as HKN [|6j, FGS [ 
[7], EKS [[8], HIJING [iS], and the most recent EPS08 [OOj. We choose the two latter 
parameterizations, which are characterized by the largest nuclear effects. Both have 
strong suppression at lower x, where we augmented HIJING with multiple scattering 
corresponding to the saturated Cronin model with C — 0.35 (GeV/c)^ as in Refs. [ 
[TTl [T^ . The EPS08 is based on the earlier EKS, however the authors included not 
only eA data, but the latest results from the RHIC dAu experiment [T^. Inclusion 
of the RHIC data results in strong suppression at small x values. In parallel a large 
anti-shadowing appears at around x « 0.2, due to the normalization. 
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Fig. 1. (Color online.) Nuclear modifications to the forward and backward direction 
based on HIJING [[9] shadowing parameterization. (See text for details.) 

Applying these shadowing parameterizations we calculated R^^uiPTi v) iii wide 
ranges of rj from backward to forward. Results are presented on Fig. [l]for HIJING. 
Forward and backward calculations are plotted with (red) dashed and (blue) solid 
respectively. We also plotted using thin dotted and dash-dotted lines on Fig. [T] 
HIJING calculations for R^A^iPTi v) without multiple scattering indicated by 'C = 
0'. The differences on Fig. [1] between the pQCD calculations with and without 
multiple scattering diverge between 77 > and 77 < 0. In the backward direction the 
lack of scattering centers in the d nucleus limits the possible collisions, thus the two 
curves merge as 77 increases. In the forward direction, without multiple scattering 
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HIJING gives a suppression at small pr related to low x values. The deviation 
is growing as 77 increases. Note the lack of experimental data in the backward 
direction. The opening with increasing \r]\ between the curves without and with 
multiple scattering suggests an 77-dependent multiple scattering similarly to Rcf. [ 
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Fig. 2. (Color online.) Nuclear modifications to the forward and backward direction 
based on EPS08 [ITO| shadowing parameterization. (See text for details.) 

On Fig. [5] we calculated the RdAuiPTiV) ^^Pplying the EPS08 shadowing pa- 
rameterizations. Solid (blue) lines represent the backward, dashed (red) lines are 
for the forward calculations. In the EKS/EPS08 frameworks the effect of multiple 
scattering (Cronin e.g. in Ref. [ 15]) is modeled with a strong anti-shadowing peak 
at around x « 0.2. At small 77 forward and backward have similar slopes, but at 
77 > 1.5 — 2.0 this trend flips over due to reaching the steep positive slope from the 
new \ow-pT RHIC dAu data. High 77 and low pr data are fitted well with EPS08, 
but at high px the anti-shadowing peak overestimates the data around midrapidity [ 

m. 



4. Analyzing Rapidity Asymmetry Data 

Fig. [3] shows the calculated rapidity asymmetry, Y^^ (pt) in the same nega- 
tive/positive 77 ranges previously presented for RdAuiPT : v) ■ Based on eq. ([4]) this 
is the ratio of the nuclear modification factors in the proper backward and forward 
rapidity ranges. Since STAR [ 2^ has published data on pseudorapidity asymmetry, 
this gives a nice opportunity to compare different models with measured values re- 
fiecting the nuclear modifications in the backward direction. EPS08, plotted with 
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(blue) dashed lines, agrees well with the data in a wide pT range. HIJING itself 
{dotted red curves) shows a similar trend, but inclusion of multiple scattering with 
C = 0.35 (GeV/c)^ turns over the slope below pt ^ 4 — 5 GeV/c at higher 77 {solid 
red lines). Similarly to the RdAuiPTiV) studies in Sec. [3l this suggests that the 
efficiency of multiple scattering should decrease going beyond midrapidity. 
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Fig. 3. (Color online.) Rapidity asymmetry calculated with HIJING [ [9] and 
EPSOS [ilOj shadowing parameterizations. (See [[2] for data and text for details.) 



5. Conclusions 

We analyzed the nuclear modification factor, R^AuiPTiV) ^^^d the (pseudo)rapidity 
asymmetry, ^^^^^(pt) in a wide rapidity range. We showed that these physical 
properties are related to each other: both reflect the geometry of the collisions and 
nuclear modifications. We found the equivalence between them in a general formula. 
This relation led us to test our calculated nuclear modification factors on backward 
hadron production, which has not been measured directly yet. 

Nuclear modifications at midrapidity follow an x scaling in a wide kinematical 
range as we presented in Refs. [[16l [17]. The pQCD improved parton model agrees 
with experimental data from CERN SPS up to RHIC energies [ [5l [12] . Several 
shadowing parameterizations and nuclear PDFs were tested to get the best x^ with 
the experimental data. We found that off-midrapidity the x-scaling is violated using 
standard shadowing parameterizations. 

To correct this problem, a rapidity dependent multiple scattering would be 
desirable, which is related to the path length (number of scattering centers) of a 
particle in the system. This would be equivalent to a strongly 6-dependent inhomo- 
geneous shadowing parameterization. 
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